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Cosmological Implications of the Very High Redshift GRB 050904 

P. A. Price, 1 L. L. Cowie, 1 T. Minezaki, 2 B. P. Schmidt, 3 A. Songaila, 1 Y. Yoshif 

ABSTRACT 

We report near simultaneous multi-color (RIYJHK) observations made with 
the MAGNUM 2m telescope of the gamma ray burst GRB 050904 detected by 
the SWIFT satellite. The spectral energy distribution shows a very large break 
between the / and J bands. Using intergalactic transmissions measured from high 
redshift quasars we show that the observations place a 95% confidence lower limit 
of z — 6.18 on the object, consistent with a later measured spectroscopic redshift 
of 6.29 obtained by Kawai et al. (2005) with the Subaru telescope. We show 
that the break strength in the R and / bands is consistent with that measured in 
the quasars. Finally we consider the implications for the star formation history 
at high redshift. 

Subject headings: cosmology: observations — gamma-rays: bursts — galaxies: 
distances and redshifts 

1. Introduction 

Identifying the epoch of reionisation remains a prime objective in modern cosmology. 
Quasars at the highest redshifts, z > 6.1 of which five are known (Fan et al. 2001, 2003, 
2004, 2005), show a complete absorption trough blueward of the redshifted Lya line from 
scattering by neutral hydrogen in the intergalactic medium (Becker et al. 2001), the so-called 
Gunn-Peterson effect (Gunn & Peterson 1965). However, because of the extreme sensitivity 
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of the attenuation to the neutral fraction of H I, it is only possible to deduce that z ~ 6 
signals the end of the epoch of reionisation (or even the end of the most recent epoch of 
reionisation; e.g., Wyithe & Loeb 2003). 

One way of robustly determining the actual epoch of reionisation is to measure the 
luminosity function of the Lya flux of high redshift galaxies, which will be attenuated by 
the damping wing of the Gunn-Peterson trough in a less severe manner than the continuum 
flux blueward of Lya. Several searches for Lya emitters (LAEs) at z ~ 6.5 have been made, 
yielding a (currently published) total of 13 that have been spectroscopically confirmed (Hu 
et al. 2002; Kodaira et al. 2003 — two; Rhoads et al. 2004; Kurk et al. 2004; Stern et 
al. 2005; Taniguchi et al. 2005 — seven new). These sources show little evolution from 
analagous populations at z ~ 5.7 (Stern et al. 2005), perhaps indicating that we have not 
yet reached the epoch of reionisation. Studies using color break samples show a drop in 
the star formation rate by about a factor of five at z ~ 6 from the peak at lower redshifts 
(Bunker et al. 2004; Bouwens et al. 2005), while the Lya surveys appear to show a flatter 
evolution (Hu & Cowie 2005). 

The color break searches are restricted to the small deep fields observed with HST, 
while, the searches for LAEs are intensive, requiring large investments of premier ground- 
based facilities, and they select only the most luminous LAEs. These surveys may therefore 
be biased towards detecting the more vigorously star-forming galaxies that produce large 
Stromgren spheres which allow the Lya flux to leak out. This means that the measurement 
of the epoch of reionisation will be biased towards higher redshifts. To address this bias, 
we require a means of homogeneously identifying more modest star-forming galaxies at high 
redshift — this is what gamma-ray bursts (GRBs) may provide. 

The luminosity function of the optical afterglows of GRBs extends as bright as an 
absolute magnitude of Mr rs —31.5 mag at 1 hour after the GRB in the rest-frame, and 
probably even brighter at earlier epochs. Because of this extreme luminosity, they can be 
detected to great distances and therefore provide an exciting way to find very high-redshift 
galaxies beyond the current upper limits of z ~ 7 and to map the star formation history at 
these extreme redshifts in a way which, if not itself unbiased, is at least independent of the 
properties of the underlying galaxies (Lamb & Reichart 2000). 

Indeed, it was widely expected that the advent of the Swift satellite would produce a 
large rate of return of very high redshift GRBs (e.g., Bromm & Loeb 2002). While the pre- 
dictions appear, in hindsight, to have been somewhat optimistic, the sensitive Swift mission 
is detecting GRBs at a higher mean redshift of z ~ 2 (Berger et al. 2005) than previ- 
ous missions such as HETE-2, BEPPO-SAX and IPN. While this may make the afterglows 
more difficult to identify, it also gave hope that SWIFT would find sources beyond the most 
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distant-known GRB (z = 4.5; Andersen et al. 2000). The burst GRB 050904 discussed in 
this paper does just that, pushing the redshift limit for GRBs beyond z — 6. Hopefully it is 
just the first of many such detections stretching to still higher redshifts. 

High redshift GRBs are easy to distinguish with coordinated optical and near infrared 
(NIR) observations, since the Gunn-Peterson effect drastically attenuates flux in the optical 
bands. The mean transmissions as a function of redshift based on quasar observations to 
z = 6.4 are tabulated in Songaila (2004), and these can be used to obtain the redshift of a 
GRB from its colors. However, Swift's Ultra- Violet and Optical Telescope (UVOT) is limited 
to observations at wavelengths bluer than about 6500A. For this reason, optical afterglows 
of GRBs at very high redshifts (z > 6) cannot be detected by the UVOT on Swift, and hence 
their properties must be characterized using ground-based NIR observations. 

GRB 050904 was triggered the Burst Alert Telescope (BAT) on Swift at 1:51 UTC on 
2005 Sep. 4 and rapidly localized (Cummings et al. 2005). Fox & Cenko (2005) undertook 
optical observations in the optical R and i' bands starting about 3.5 hours after the GRB, but 
did not identify any afterglow candidate to reasonable limiting magnitudes. The subsequent 
identification of a bright afterglow in the NIR J band led to the interpretation that this was 
a GRB at very high redshift (Haislip et al. 2005). 

Observations of the afterglow with the MAGNUM telescope began about 12 hours after 
the GRB (§2). These observations allow us to generate a spectral energy distribution for 
the source at that time which yields a strong lower limit on the redshift of z = 6.18. This is 
consistent with the photometric redshifts reported by Haislip et al. (2005) and Tagliaferri 
et al. (2005), and the spectroscopic redshift of z — 6.295 ± 0.002 measured by Kawai et 
al. (2005). Furthermore, the limits on the break strengths R — J and I — J are consistent 
with the object lying at the spectrosopic redshift (§3). We describe these observations in the 
present paper and briefly speculate on the implications for the star formation history of the 
universe (§4). 

2. MAGNUM observations of GRB050904 

MAGNUM (Multicolor Active Galactic NUclei Monitoring) is a 2 meter telescope on 
Haleakala built by the Research Center for the Early Universe (RESCEU) at the University of 
Tokyo and used to study AGN variability (Yoshii 2002). In order to optimise the efficiency 
of the monitoring observations, the telescope is operated in a robotic mode using queue 
scheduling. GRB observations can be performed as soon as a notification is received by 
inserting the target and overriding the queue. 
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MAGNUM's principal instrumentation is the Multicolor Imaging Photometer (MIP; 
Kobayashi et al. 1998), a dual-beam optical/NIR camera which covers a 1.5 arcmin square 
field in the UBVRIY JHK S KL' bands (though U observations are difficult, and V infea- 
sible). The instrument, mounted at the Cassegrain focus, uses an internal beam splitter to 
send the short wavelength light to a 1024 pixel square thinned CCD (though the entire CCD 
is not illuminated), and the long wavelengths to a 256 pixel square InSb array. 

Because the limited field of view of the MIP makes it impractical to observe the 4 arcmin 
localisations from Swift's Burst Alert Telescope (BAT), we target afterglows discovered by 
the X-Ray Telescope (XRT) or ground-based follow-up observations in order to characterise 
their spectral flux distribution, and so attempt to determine a photometric redshift from the 
Lya absorption at high redshift. We use a pre-planned sequence of four observations, each 
consisting of nine individual minute-long exposures with a box dither pattern of 10 arcsec 
step: RIRI in the optical and HKJY in the NIR. The final images provide an accurate and 
nearly simultaneous spectral energy distribution for the object. The entire sequence takes 
about 84 minutes, with an additional 5 minutes before commencing the GRB observations 
in order to correct the pointing and focus after the slew. 

We observed GRB 050904 in this manner between 13:57 and 15:07 UTC on 2005 Sep. 
4, or about 0.51 days after the GRB. The final combined images are shown in Figure 1. 
Because the only 2MASS source in the field is a galaxy, the images could not be immediately 
photometrically calibrated. We obtained observations on subsequent nights to obtain UK 
calibrations of two stars in the field. The H band was calibrated from the 2MASS galaxy 
by using a large aperture to measure the entire flux of the galaxy, and applying a measured 
aperture correction to obtain the flux of the afterglow. The R band was calibrated using 
SDSS magnitudes for a star in the field, and applying the appropriate transformation (Smith 
et al. 2002). The uncertainties in these two non-ideal calibration methodologies do not 
strongly affect our results, since the R and H measurements are less important than the 
UK measurements. The fluxes in /xJy with 1 a errors are summarized in Table 1. These 
fluxes are measured in matched apertures corrected for the image quality, with radii ranging 
from 1.4 to 1.7 arcsec. The Y band observations are not particularly sensitive, so are not 
included in our analysis. Since no source is apparent in the RI bands, we measured the flux 
in the image at the position of the afterglow. Errors were taken as the r.m.s. from zero of 
the distribution of flux in apertures randomly distributed over empty background regions. 

The difference between the J-band magnitude reported here and that obtained by 
UKIRT at a slightly earlier epoch, combined with the unusual behaviour in the Z band 
at about this epoch reported by Haislip et al. (2005) likely indicates that the source was 
somewhat variable during the period of these observations. Such short timescale variability 




Fig. 1. — Images of GRB 050904 (circled object) in the RIYJHK bands respectively. Each 
image is 1.2 arcmin on a side. The detector response is poor at Y and the failure to detect 
the afterglow at this wavelength is not significant. 
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has been detected in other GRB optical afterglows, with GRB 030329 a notable example (Lip- 
kin et al. 2004). This serves to demonstrate the need for simultaneous or near-simultaneous 
multicolor imaging, such as MAGNUM provides. 

3. Redshift limits and the break strength 

We show the measured spectral flux distribution (SFD) of the source in Figure 2. There 
is a substantial break at the / band which places a tight lower limit on the redshift of the 
source. Because the I band data is consistent with a null detection we cannot place any 
useful upper limit (better than z <^ 8) on the redshift of the source. 

In order to obtain the redshift estimate, we fitted a power-law spectrum to the JHK 
data, obtaining a spectral slope of f3 = 0.3 ± 0.6, where f v oc u"^; though the error bar is 
large, such a shallow spectral slope is likely a product of intrinsic variability in the source 
during the J-band observation, as discussed earlier. We then modulated this spectrum with 
the Lya and Ly/3 transmissions of the intergalactic medium measured by Songaila (2004) 
in high redshift quasars. The resulting spectrum is shown for z = 6.18 in Figure 2, where 
we also show the positions of the / filter. In order to reduce the / band flux to observed 
value we require z > 6.18 at the 2 a level and z > 6.37 at the 1 a level. The results are 
extremely sensitive to the adopted redshift. They are consistent with the redshift range of 
z = 6.30 ± 0.07 reported by Tagliaferri et al. (2005) based on similar observations with the 
VLT and also with the spectroscopic redshift of z = 6.295 ± 0.002 found by Kawai et al. 
(2005). 

In Figure 3 we show the break strength between R and J, and between / and J, directly 
compared with the measured values in individual quasars at these redshifts. The quasar 
values are measured by comparing a power law fit to the continuum in line free regions of 
the quasar to the directly measured flux in the / band (Songaila 2004). The R band provides 
a weaker constraint than the / band but would still place the GRB at z > 6.1. 

4. Discussion 

Observations of high redshift GRBs can be used to infer the cosmological star formation 
history either through observations of the host galaxy or by translating the GRB rate to a 
star formation rate (Totani 1997). The most direct method is to assume that the rate of 
GRBs as a function of redshift is proportional to the rate of formation of high-mass stars and 
hence (assuming that the mass function is invariant) to the total star formation rate. In order 
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Fig. 2. — The spectral flux distribution of GRB 050904 from the MAGNUM observations 
are shown by the solid squares with 1 a error bars. The dotted line shows the position of the 

1 filter. The dashed line shows the expected SFD of the object at a redshift of z = 6.18 (our 

2 a lower limit on the redshift) based on a power law fit to the longer (JHK) wavelengths, 
modulated by the measured transmission of the intergalactic medium below the redshifted 
Lya position computed using the measurements of Songaila (2004). 
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Fig. 3. — 2 a lower limits on the I — J and R— J break measured in the GRB (solid diamonds) 
are compared with directly measured values in high redshift quasars (solid squares). 
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to make such an interpretation, we need to calibrate the GRB rate versus the star formation 
history at lower redshifts, and it will always be subject to the assumption that the fraction 
of massive stars forming GRBs and the initial mass function of the stars remain invariant at 
the higher redshifts (e.g., Porciani & Madau 2001; Lamb & Reichart 2000; Bromm & Loeb 
2002; Natarajan et al. 2005). 

Prior to GRB 050904 twenty spectroscopic redshifts had been measured for Swift GRBs 
(excluding short/hard GRBs), of which fourteen lie beyond z = 1. Thirteen of these fourteen 
would still have been detected above the Swift BAT threshold of 0.2 ph/cm 2 /s if they had 
lain at z — 6.29. This indicates that the selection effects, at least in the gamma-rays, are 
not strong. Here we make the simple assumption that the thirteen GRBs summarized in 
Table 2 which would have fluxes about 0.2 ph/cm 2 /s at z = 6.29 represent the low redshift 
counterparts of GRB 050904 detected over the same period of Swift observations, which 
indicates that the efficiency of detecting GRBs at z = 6.29 is near unity, compared to GRBs 
at lower redshifts. For each GRB we give the redshift, the observed peak flux and the value 
of the peak flux if the source had been at z = 6.29. 

We point out that this is, of course, a simplified analysis and that there are more 
selection effects than just the detection of the GRB itself. In particular, the success rate 
of detecting the optical/NIR afterglow will be a function of redshift, as will the fraction of 
afterglows for which it is possible to measure the redshift, either from absorption or emission 
lines. This is a complicated endeavour requiring detailed Monte Carlo simulations of GRB 
afterglow searches to determine which GRBs with and without redshifts might plausibly 
have been identified if followed up in the same way as GRB 050904, and is beyond this 
simplified analysis. However, such an analysis is not justified at the present time, given the 
significant small number uncertainties in the data and the desire for a larger sample of Swift 
events to use in bootstrapping the completeness estimate. Nevertheless, our measurement 
will provide, at the least, a lower limit to the star formation density, since we know that we 
are missing some GRBs. 

Because of these complicated selection effects, which have discouraged all but the most 
bold (e.g., Blain & Natarajan 2000) from using the actual redshifts measured from the 
optical afterglow or host galaxy in determining the GRB rate as a function of redshift, most 
other attempts have relied upon assuming that GRBs are standard candles (e.g., Totani 
1999), or using empirical "pseudo-redshift" indicators from the GRB itself (e.g., Firmani 
et al. 2004; Lloyd-Ronning, Fryer & Ramirez-Ruiz 2002) such as the luminosity-variability 
relation (Reichart et al. 2001). These methods buy a large sample size and better understood 
selection effects at the cost of uncertain (by at least a factor of two!) redshifts. This results 
in a need to deconvolve the resulting rate distribution using the redshift errors, making it 
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difficult to determine the true behaviour of the GRB rate at high redshift . 

The relation between the star formation history and the GRB rate has been derived by 
a number of authors (e.g. Porciani & Madau 2001; Lamb & Reichart 2000; Bromm & Loeb 
2002; Natarajan et al. 2005). At z > 1 where the effects of the cosmological constant are 
negligible, the number of gamma ray bursts per unit redshift, dN/dz, in a given observing 
time interval is simply related to the star formation rate per unit comoving volume, ip, by 
the function: 

dN/dz — Aip [(1 + z)~ 2 - 5 - 2(1 + z)~ 3 + (1 + ^)~ 3 - 5 ] (1) 

where the normalizing factor, A, is assumed to be independent of z. 

In Figure 4 we compare the shape of the star formation rate determined from the 
Swift GRBs with color-selected galaxy determinations of the star formation rate over the 
same redshift range taken from the paper of Bunker et al. (2004). We have set A = 
0.005 M Q Mpc~ 3 to match the observations at z — 3. Within the wide uncertainties left by 
the small number of statistics, the current values cannot differentiate between the slow decline 
seen in the color selected galaxies and a flat star formation rate with redshift. However, it is 
clear that as the sample size increases we should be able to make a valuable comparison. The 
GRB determinations are more powerful in some ways, since they relate to all star formation 
including those in lower luminosity galaxies than can be directly detected at the present 
time. This holds promise that future identifications of z > 6 GRBs will enable a complete 
measurement of the star formation rate density at very high redshift. 



5. Summary 

In the present paper we have described the observations of GRB 050904 with the MAG- 
NUM telescope. These observations place a strong lower limit on the redshift of z = 6.18 
consistent with the spectroscopic redshift of z = 6.29 measured by Kawai et al. (2005). 

The most immediate result is that GRBs exist at z > 6, and that they can be identified 
using a simple set of near-simultaneous optical and NIR observations. This presents the 
prospect of using the afterglows of high-redshift GRBs not only as lighthouses to illuminate 
the high-redshift Universe (as is done for quasars today), but also as signposts to alert 
observers to the presence of the host galaxy, allowing deep follow-up observations to measure 
the Lya flux. With the discovery of more high-redshift GRBs, it should be possible to form 
a useful sample for nailing down the epoch of reionisation. 

We also gave a simple discussion of the star formation rate history from z = 1 — 7 based 
on the current Swift GRB observations showing that within the still broad uncertainties, the 
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Table 1. Measured fluxes from MAGNUM observations in multiple bands. 



Band 


Flux (/iJy) 


R 


-1.04 ±0.75 


I 


0.4 ± 1.2 


J 


47.9 ±6.2 


H 


30.0 ±8.9 


K 


41 ± 14 



Note. - The observations were made at a mean epoch of 2005 Sep. 4.60, and may be 
treated as simultaneous — any correction to the fluxes for the decay of the afterglow over the 
course of the observations would be smaller than the measurement errors (0.15 mag, using a 
temporal decay index of 1.2; Haislip et al. 2005). These measurements are not corrected for 
the relatively small foreground Milky Way extinction. The R and I band measurements are 
consistent with no detections to 3a upper limits on the fluxes of 2.3 and 3.6 /iJy, respectively. 
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Fig. 4. — The star formation rate as a function of redshift in units of solar masses per 
Mpc 3 , taken from the compilation of Bunker et al. (2004) is shown by diamonds. The solid 
boxes with 1 a error bars show the corresponding determinations from the Swift GRBs for 
a normalizing constant of 0.0033 in the same units. 
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Table 2. Lower redshift Swift GRBs counterparts to GRB050904 



GRB 


Redshift 


Peak flux 


r 

7.2=6.29 


References 


051111 


1.549 


2.50 


0.25 


(1), (2) 


051109A 


2.346 


3.70 


0.77 


(3), (4) 


050922C 


2.198 


7.36 


1.37 


(5), (6) 


050908 


3.350 


0.70 


0.27 


(7), (8) 


050820A 


2.612 


2.50 


0.63 


(9), (10) 


050802 


1.710 


2.65 


0.32 


(11), (12) 


050730 


3.969 


0.57 


0.28 


(13), (14) 


050603 


2.821 


27.6 


7.89 


(15), (16) 


050505 


4.270 


1.81 


1.01 


(17), (18) 


050401 


2.900 


12.6 


3.77 


(19), (20) 


050319 


3.240 


1.45 


0.52 


(21), (22) 


050318 


1.440 


3.20 


0.28 


(23), (24) 


050315 


1.949 


1.98 


0.30 


(25), (26) 


050126 


1.290 


0.70 


0.05 


(27), (28) 



Note. - The sample is limited to Swift GRBs at 2 > 1. The peak fluxes are in 
ph/cm 2 /s, measured in the 15-150 keV band; they are taken from the Swift archive 1 . All 
but GRB 050126 would have been detected by Swift if placed at z = 6.29. References: (1) 
Hill et al. (2005); (2) Krimm et al. (2005); (3) Quimby et al. (2005); (4) Fenimore et al. 
(2005); (5) Jakobsson et al. (2005); (6) Krimm et al. (2005); (7) Fugazza et al. (2005); (8) 
Sato et al. (2005); (9) Prochaska et al. (2005); (10) Cummings et al. (2005); (11) Fynbo 
et al. (2005); (12) Palmer et al. (2005); (13) Chen et al. (2005); (14) Markwardt et al. 
(2005); (15) Berger & Becker (2005); (16) Fenimore et al. (2005); (17) Berger et al. (2005); 
(18) Hullinger et al. (2005); (19) Fynbo et al. (2005); (20) Sakamoto et al. (2005); (21) 
Fynbo et al. (2005); (22) Krimm et al. (2005); (23) Berger & Mulchaey (2005); (24) Krimm 
et al. (2005); (25) Kelson & Berger (2005); (26) Krimm et al. (2005); (27) Berger, Cenko 
& Kulkarni (2005); (28) Sato et al. (2005). 
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observations point to a flat or a slowly declining star formation rate consistent with color 
selected galaxy observations. 

GRB 050904 is an exciting precursor to further high redshift GRBs which should allow 
us to refine the star formation analysis and to study the properties of the intergalactic 
medium at these redshifts through color break and spectroscopic techniques. 

We thank Elizabeth Stanway for providing the data for Figure 4 in tabular form. This 
work was supported by a Swift Guest Investigator grant (NNG05GF40G) and a Grant in Aid 
of Center of Exellence Research (07CE2002) of the Ministry of Education, Science, Culture 
and Sports of Japan. We thank the anonymous referee for helpful suggestions that improved 
this paper. 
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